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Improved method and device for local spectral analysis of an ultrasonic signal 
DESCRIPTION 

Technical field 

The present invention relates to a method to process ultrasonic signals, 
5 especially in the field of non-invasive clinical tests based on the ultrasonic 
method. The Invention also relates to an ultrasound device for the 
implementation of a processing method of the ultrasonic signals. 
State of the art 

Ultrasound scanning is a non-invasive investigation method for 

10 examining internal organs of the human and animal body and also for non- 
invasive examination of other types of structures, including non-living and 
inorganic structures. As is known, it is based on the production of a series of 
ultrasonic signals (forming an acoustic beam) which are directed into the 
patient's body and which excite an ultrasonic response signal, caused by the 

15 back reflection and/or scattering of sound waves by the diffusers and/or 
reflectors present within the region under observation. The probe for 
transmitting and receiving the ultrasonic signal sends the response signal to a 
receiving and processing system, after which the signal is transformed into a 
two-dimensional ultrasound image. 

20 Commercial ultrasound scanners supply images relating to the 

amplitude of the ultrasonic signal back-propagated by the various tissues or 
interfaces struck by the acoustic beam. 

The various tissues or the various interfaces encountered are detected 
on the basis of their echogenicity, determined by their acoustic impedance, 

25 which in turn is the product of the density and the acoustic velocity. 

The purpose of conventional ultrasound scanners is to reconstruct and 
maintain the spatial order of the tissue under examination, by means of the 
amplitudes detected in the time domain, which in turn provides the spatial 
location in depth. The conditioning performed on the signal is aimed at 

30 maximizing the signal to noise ratio and keeping the sensitivity of the 
investigation constant when the depth varies. 

When the acoustic wave is propagated in a medium and encounters a 
non-homogeneous region, some of the energy is back-propagated towards 
the transducer and some is propagated into the tissue. 
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The back-propagated acoustic signal consists of the temporal 
overlapping of the contributions from each reflector and/or diffuser; the result 
is a signal having an extremely complex shape, due to the random nature of 
the distribution and the shapes and mechanical properties of the diffusers 
5 and reflectors. The acoustic signal back-propagated by the medium under 
examination is converted by the transducer into an electric signal, called radio 
frequency signal. 

Figure 1 shows a portion of a radio frequency signal forming part of an 
ultrasound trace relating to biological tissue. More specifically, Figure 1 

10 shows a portion of a trace of the return ultrasonic signal from a portion of an 
organ subjected to ultrasound examination. The response time, that is, the 
depth of the location of the diffuser or reflector which produces the return 
signal captured by the trasducer, is shown on the abscissa. The amplitude of 
the signal is shown on the ordinate. The ultrasonic image of the organ under 

1 5 examination is reconstructed from a set of traces or tracks of this type, which 
form a frame of signals. These traces or tracks are referred to sometimes 
also as scan lines. 

A first very simple processing operation carried out on the radio 
frequency signal is extraction of the signal amplitude by means of an 

20 envelope extraction operation. In Figure 2, the letter I indicates the result of 
the envelope operation carried out on the radio frequency signal shown in 
. Figure 1. It will be noted that the envelope signal follows the positive crests of 
the radio frequency signal and therefore provides an estimate of the signal 
amplitude. 

25 The envelope operation is currently used in all ultrasonic imaging 

techniques and forms the basis of the A-mode presentation (where A stands 
for amplitude) of the individual trace and B-mode presentation (where B 
stands for brightness) for the resulting images, the amplitude of each 
individual A-mode trace being encoded on a grey or color scale. 

30 In ultrasound scanners of the present generation, the envelope of the 

ultrasonic signal is used, that is, the detected signal is represented in the A- 
mode, B-mode and M-mode traces. 

With this solution, there is a loss of the information contained in the 
radio frequency signal, that is, in the form produced by the receiving 
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transducer, consisting of a signal having a carrier frequency approximately 
equal to that of the center frequency of the band of said receiving transducer, 
with a frequency modulation (asymmetrical spectrum) with phase steps 
(frequency peaks) and possibly with a morphology due to non-linear distortion 
5 when the medium shows non-linear effects (production of harmonic spectra). 

Figure 3 shows the envelope signal of the radio frequency trace in 
Figure 1, in a conventional A-mode representation. It will be noted that 
information on the frequency and phase of the original signal has been 
irretrievably lost. 

10 In particular, if the radio-frequency is eliminated by an envelope 

extraction operation, there is a loss of: 

• information obtained from non-linear interactions between ultrasound 
and biological tissues (local acoustic impedance which depends on the 
density of the medium, which in turn depends on the intensity of the 

15 incident ultrasonic wave, re-emission on the harmonics and 

subharmonics by bubbles present in the contrast medium, etc.); 

• information obtained from linear interactions, that is, modifications of the 
signal spectrum due to selective absorption, especially where this is 
depth-dependent, and the scalloping effect (indentation of the shape of 

20 the amplitude spectrum present in signals which have been 

backscattered by discontinuities, including those of the point type, 
distributed periodically in space). 

Conventional ultrasound scanning systems used for medical diagnosis 
are therefore based on the observation of an A-mode trace or a B-mode 
25 image, which show only the amplitudes of the echo signal relating to the 
various tissue structures encountered. 

Consequently there are situations in Which the conventional methods 
are found to be limited since, owing to low sensitivity or low specificity, it is 
not possible to obtain a reliable diagnosis of pathology. However, a radio 
30 frequency signal contains not only the information relating to amplitude, but 
also the frequency and phase information (related to the time delays with 
which the echoes due to individual diffusers reach the transducer). 

The latest signal processing methods tend to extract all the parameters 
characterizing the signal, in particular what is known as the spectral content. 
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Spectral analysis applied directly to the radio frequency signal can extract 
information from all the characteristics of the echo signal as it is returned 
from the transducer. 

A spectral analysis algorithm for "reading" the signal in the frequency 
5 domain makes it possible to provide characteristic values for differentiation of 
the various types of signal and consequently of the biological structures with 
which they interact. 

Diagram (A) in Figure 4 again shows the radio frequency signal in 
Figure 1, while the diagrams (B) and (C) show the spectral content of the 
10 signal for two portions of the trace in Figure 4A, to show how the structural 
differences in the tissues result in different types of ultrasonic impulses and 
consequently distinctive frequency variations in the respective spectrum. The 
spectral contents of the two portions of signal shown in Figures 4B and 4C 
extend from 0 to 12 MHz. It will be noted that both spectra have 
1 5 approximately the same central frequency (the frequency at which the peak is 
found) in the vicinity of 4 MHz (the signal relates to scanning with an 
ultrasonic probe having a nominal frequency of 5 MHz), but have a different 
frequency variation, since the first graph has contributions (above 5 MHz), 
which are present to a lesser extent in the other portion of the signal. 
20 In practice, the ultrasonic signals backscattered from the 

microstructures, due to the differences in the propagation path, interfere on 
the face of the receiving transducer in a destructive or constructive way, 
resulting in the "speckle" effect in the images. The difference in the amplitude 
of the echo signal, that is, the echogenicity, depending on the type of 
25 interference, is also influenced by the reflectivity of the tissues which in turn is 
determined by the mechanical nature, density and elasticity, of the 
components of the microstructure; for example, a high content of collagen 
increases local mechanical impedance and thereby echogenicity. Therefore, 
the following contribute to the appearance of the classic ultrasonic image: 
30 spatial agglomeration of the microstructures, which determines the speckle 
characteristics, and their mechanical nature which determines echogenicity. 
The mean amplitude variation in space therefore seems not to provide 
univocal characterization of the tissues as it is influenced by various 
morphological factors which are formed randomly. 



-5- 



FC1/1I2U04/UUU548 



If information can be extracted on the speckle in terms of spatial 
periodicity (spatial frequency) then it is possible also to acquire information 
related to any spatial periodicity of the agglomerations under examination. 
The radio frequency ultrasonic signal is influenced by the information 
5 produced by the speckle of the images, in terms of particular amplitude and 
phase modulations of the ultrasonic signal. 

The information present in the time domain are preserved in the 
frequency domain and the influence of the speckle Is translated into the 
scalloping effect that consists in a variation of amplitude according to the 
10 graphical representation indicated in Figure 4D. In fact, for example in the 
case of N diffusers or reflectors (the distinction, which is not unequivocal, is 
made in relation to the wavelength of the ultrasonic signal) struck by the 
ultrasonic beam, the spectral density of power, calculated according to the 
Fourier transform, of the echo signal received |x(<y)| 2 may be represented as: 

15 \x(pf -Istorj + 2% ± Vj cos«,(2(/, -/,)/,)] 

where: 

W^l : spectral density of power of the signal transmitted 
ri: reflection coefficient of the i-th reflector 
li: distance between the first and the i-th reflector 
20 c: velocity of the ultrasounds. 

In the case of only two reflectors separated by a distance equal to 1, the 
expression takes a simpler form: 

\X(4 2 = \s(4 2 [ n 2 + r 2 2 +2 m cos(2^/c)] 
T = = 7 ' S - he distance in time betwee " the peaks of the signals 
25 back-propagated by the two reflectors. 

It can be seen from the schematic formulae that the spectrum is 
weighted by a cosine with period proportional to the distance of the reflectors; 
this determines the minimum and maximum values of the spectral power 
density. This effect is called scalloping, and is represented schematically in 
30 Figure 4D for two reflectors, positioned at variable distances (indicated with 
SO, S1 , S2) from each other in the direction of propagation of the ultrasonic 
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beam. 

In order to distinguish two otherwise indistinguishable reflectors, that is, 
at distance smaller than the duration of the impulse transmitted, it is possible 
to evaluate the distance between minimums in the spectral domain. In the 
5 case of several reflectors the spectrum of the initial signal is weighted by 
overlapped oscillations having different periods. The formulae also include 
the coefficients of reflection of the reflectors which determine the amplitude of 
the backscattered signal. 

From the formula it can be seen that the information present in the time 
10 domain influence the shape of the spectrum, so that an accurate description 
of the variation in distribution of the amplitudes in this domain should supply 
the elements for differentiation of the biological structures. 

These exemplificative considerations show how spectral analysis on the 
radio-frequency signal is a powerful instrument for the characterization of 
1 5 tissues by means of the relative ultrasonic signal. 

In recent years the development of new signal processing algorithms 
has made it possible to obtain further information on tissues, by investigating 
interaction with the ultrasonic beam in terms of spectral variations induced in 
the reflected signal. In fact, it has been found that pathological structures can 
20 be identified according to the differences in their selective frequency filtering 
action. 

However, the methods developed up to the present time have not 
proved to be effective, either because they are unsuitable to characterize 
tissue reliably, or because of the excessive computing capacity required for 
25 the algorithms used, which cannot be implemented in a way that permits real 
time display of the results and would therefore not allow clinical applications 
where the speed of data output is a basic requirement 

EP-A-1.341.003 describes a new method for spectral analysis of an 
ultrasonic signal, which overcomes the limits and drawbacks of prior art and 
30 allows extraction of the information producing the speckle of the images. The 
content of this previous publication is fully incorporated in the present 
description and forms an integral part thereof. 
Objects and Summary of the Invention 

Starting from the method described in EP-A-1.341.003, the object of the 
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present invention is to provide an improvement and extension to the 
possibilities offered by this type of spectral analysis, to obtain greater and 
more precise information on the structure under examination. 

Essentially, starting out from the same premises on which the method 
5 described in EP-A-1 .341.003 is based, the invention described herein allows 
the description of the spectrum to be improved by means of: a suitable choice 
of the mother wavelet and/or the choice to decompose the entire spectral 
band of the signals into a variable number of sub-bands of arbitrary width and 
position, so as to extract from the sample under investigation the greatest 
10 possible amount of information and/or a more effective use of the coefficients 
characteristic in the description of the local spectrum. 

In practice, according to a first aspect, the invention is based on the 
idea of combining local estimators, hereinafter also called local spectral 
parameters, obtainable from the ultrasonic signal by means of a procedure 
15 analogous to the one described in EP-A-1 .341 .003, using information relative 
to the shape of the distribution histograms in various classes of values of said 
local estimators. In practice, from a frame of radio frequency signals collected 
from a structure under examination, one or more frames or matrices of local 
spectral parameters or local estimators are obtained. Subsequently, these 
20 estimators are processed statistically, by moving a window with pre- 
established dimensions along the matrix of the local estimators to determine 
the distribution histogram of the local estimators in various classes of values. 
A parameter relating to their shape, for example the standard deviation, is 
obtained from these histograms. The local estimators are then combined, for 
25 example multiplied, divided or processed in another way together with the 
shape parameters of the distribution histogram. A new matrix of local 
weighted estimators is thus obtained, to allow value intervals of the 
estimators to be associated with homogeneous regions of the investigated 
sample. These regions will then be identified by histological analysis. 
30 Representation of these regions of the original ultrasound frame takes place 
by means of a color code overlaid on the conventional B-mode image. 

The method of the present invention may be abbreviated to the RULES 
method (Radiofrequency Ultrasonic Local Estimator) as it is based on the use 
of local estimators which contain spectral information regarding the ultrasonic 
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radio frequency signal returned from the structure examined. 

As shall be apparent from the text hereunder, the local estimators are 
represented by the coefficients of a suitable degree, for example fourth or 
fifth degree, interpolating polynomial, which approximates the coefficients 
5 obtained for the various frequency bands into which the digitized and 
sampled radio frequency signal was previously subdivided. Alternatively, the 
iocal estimators may consist in combinations of coefficients of various 
degrees of the interpolating polynomial. The combination of coefficients of 
various degrees can take place before or after having weighted each 
1 0 coefficient with the shape factor of the distribution histogram. 

According to a different aspect, the invention consists in processing the 
coefficients of the interpolating polynomial to obtain local estimators from a 
combination of two or more coefficients of different degree of said 
polynomial, without weighting the coefficients by the shape parameter of the 
15 distribution histogram. Also in this case greater information is obtained with 
respect to prior art methods, as several combined coefficients of the 
interpolating polynomial are utilized instead of a single coefficient. 

Subdivision into bands, that is, decomposition into bands or sub-bands 
can take place with various criteria. Moreover, also the shape parameter of 
20 the distribution histogram of the local estimators can be chosen in various 
ways. In particular, it is possible to; 

1 . use different Mother Wavelets to adapt to the typical speckle of 
the investigated tissue; 

2. use different levels of decomposition both to obtain different 
25 spatial resolutions and to use various subdivisions of the overall spectrum 

available in a combined way; 

3. use combinations of the coefficients of the interpolating 
polynomial "weighted" with parameters related to the nature of their statistical 
distribution by means of a targeted study of the histograms through specific 

30 parameters such as: 

most populated classes 
standard deviation 
kurtosis 

higher degree moments 
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symmetry index 

The combination of suitably weighted coefficients produces, for each 
window (the dimensions of which can be set) of the processed ultrasound 
image, a new matrix consisting of these combinations. A new statistic is 
5 performed on it and the classes associated with the various structures to be 
characterized are chosen. To increase the differentiation capacity of the 
tissues, that is, identification of the "typical" classes, combined analysis of the 
classes is once again used. 

Particularly advantageous embodiments of the method according to the 
10 invention are indicated in the attached claims and shall be described in 
greater detail below with reference to some examples of embodiments and 
experimental tests performed. 

According to a different aspect, instead of performing a combination or 
weighting operation of the local estimators and the shape parameters of the 
15 distribution histograms of said estimators, according to the invention two or 
more local estimators can be combined and this combination can be used as 
a combined local estimator. For example, several matrices of local estimators 
can be produced, in which the individual coefficients are constituted, in each 
matrix, by coefficients of various degrees of the interpolating polynomial. 
20 Subsequently, coefficients in homologous positions of these matrices can be 
combined, for example the coefficients of all or some of said matrices can be 
added or multiplied together. In conclusion, in this way a matrix of coefficients 
is obtained, each of which is composed of a function of the local estimators 
contained in the previously obtained matrices. 
25 Brief description of the drawings 

The invention will now be more clearly understood from the description 
and the attached drawing, which shows a non-limiting practical embodiment 
of the invention. More specifically: 

Figures 1 to 4A-4D, already described above, show an ultrasound trace 
30 and a representation of the spectral content in two portions of said trace; 

Figure 5 shows a simplified block diagram of an ultrasound scanning 
apparatus; 

Figure 6A shows an ultrasonic excitation signal sent by the transducer 
into the organ to be examined; 
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Figures 6B to 6D respectively show the RF response signal, the 
Daubechies 16 wavelet used for filtering the digitized RF signal and the 
coefficients of the DWPT obtained by decomposition into eight bands; 

Figures 7A-7I schematically show the processing sequence on the 
5 ultrasonic signal performed by the method according to the present invention; 

Figures 7J, 7K schematically show modified embodiments of the 
processing procedure according to the inventions; 

Figures 8A and 8B show the B-mode ultrasonic images respectively of 
the carotid artery and prostate, on which the results shown in Figures 9A to 
10 1 7B were obtained ; 

Figures 9A-9C show the variation of three fourth degree interpolating 
polynomials in a particular application to the examination of a carotid artery, 
relating to the portion of lipid plaque, blood and calcification; 

Figure 10 shows the distribution histograms of the various orders of 
15 coefficients of the fourth degree interpolating polynomial relating to the 
portion of lipid plaque of the carotid artery; 

Figure 1 1 shows the distribution histograms of the coefficients of the 
interpolating polynomial in the blood; 

Figure 12 shows the distribution histograms of the coefficients of the 
20 interpolating polynomial in the calcification region; 

Figures 13A-13C show the variation of the fifth order interpolating 
polynomials obtained by the same application as in Figures 9A-9C; 

Figures 14, 15 and 16 show, analogous to Figures 10, 11 and 12, the 
distribution histograms of the various coefficients of the fifth degree 
25 interpolating polynomial in the various regions of lipid plaque, blood and 
calcification; 

Figures 17A and 17B show the distribution histograms of the 
coefficients of the fourth order interpolating polynomial applied to the analysis 
of a prostate gland; 

30 Figure 18A shows a block diagram of a cascade of filters for 

decomposition of the digitized and sampled radio frequency signal into six 
bands or sub-bands by means of a discrete wavelet, in bands with variable 
amplitude, with increased subdivision in the low frequency regions; 
Figure 18B schematically shows the six decomposition bands; 
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Figure 18C shows the variation of the fourth order interpolating 
polynomial obtained by applying the decomposition in Figures 18A, 18B to an 
ultrasonic response signal of a carotid artery; and 

Figures 19, 20 and 21 each show a B-mode ultrasound image of a 
5 portion of a carotid artery affected by lipid plaque and calcification 
phenomena, and the variation of the interpolating polynomial and the 
distribution histograms of the coefficients of the various degrees of the 
interpolating polynomial, in the three regions affected respectively by lipid 
plaque, blood (free vessel) and calcification; 
10 . Figure 22 shows a structure of the decomposition tree obtained by a 
DWPT transform; and 

Figures 23A and 23B show explanatory graphs of the classification 
procedure of the local estimators. 

Detailed Description of the Preferred Embodiment of the Invention 

15 Figure 5 shows a highly schematic block diagram of an ultrasound 

scanning apparatus. The number 1 indicates a transducer of the scanning 
head, indicated generically with 2, which emits the ultrasonic excitation signal 
and receives the ultrasonic response signal. The analog output 3 of the 
scanning head is connected to an acquisition card 5. The latter produces a 

20 sampled and digitized radio frequency signal, which forms the input for a 
bank of integrated FIR filters, indicated generically with 7. The bank of filters 
7 can consist of any commercially available integrated device, for example 
the 43168 circuit produced by Harris, the GF191 circuit produced by 
Gennum, the GEC 16256 circuit produced by Plessey, or the GC2011 circuit 

25 produced by Graychip. The FIR filter 7 is associated with a digital signal 
processor (DSP) indicated with 9, which acts as a controller for the whole 
processing unit. The microcontroller 9 can be, for example, a TMS 32031 
DSP produced by Texas (USA) or another equivalent. 

The number 1 1 indicates a memory associated with the filter 7 and with 

30 the controller 9. The memory 1 1 is connected to an analog/digital converter 
13 and to the digital stages of the scanning head 2. 

The circuit described here is known and does not require more detailed 
description. Essentially, the apparatus operates as follows. The transducer 1 
emits a series of ultrasonic excitation signals having, for example, the form 
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shown in Fig. 6A, into the body containing the organ to be examined, for 
example the prostate, the liver or another part. The tissue struck by the 
ultrasonic waves of the signal emitted by the transducer 1 returns a scattered 
or reflected signal which has the form of the signal in Fig. 1 (radio frequency 
5 signal). 

The radio frequency signal is taken from the receiving circuit of the 
ultrasound scanner immediately after the time gain compensation (TGC) unit 
and is then digitized at a suitable sampling frequency, for example 40 MHz, 
and with a resolution of 12 bits (or other suitable resolution), by the 

10 acquisition card 5. In practice, acquisition is carried out for individual traces 
until a frame or a portion of a frame is completed. The signal samples 
acquired for each individual trace are stored in a memory bank of the 
acquisition card until a complete frame or a portion of a frame has been 
acquired; it is only at this point that the stored data are sent to the processor 

15 card comprising the DSP 9 and the FIR filter 7 which executes the operations 
to be described in detail below, according to the method of the present 
invention. While the data is being transferred to the processor card, the 
acquisition card starts to acquire the signal samples of the traces for the next 
frame. The processing time of the processor card 12 is short enough to allow 

20 processing to be completed within the acquisition time of the sampled signal 
for the next frame. Thus real time operation is achieved. 

the processor card 12 carries out processing of the signal, particularly 
the time-frequency transform using the integrated filters 7 and the other 
operations to be described in greater detail below, to obtain characterization 

25 of the signal spectrum. 

Figures 7A-7E summarize, in the form of graphs, the operations carried 
out on the input signal to obtain the representation of the matrix containing 
the information relating to the suitable combination of local estimators 
determined using the method according to the present invention. For a 

30 detailed description of the operations performed in the various processing 
phases in a specific example, reference can be made to the description in 
EP-A- 1.34 1.003, where not described in detail herein. 

Processing of the sampled and digitized signal acquired by the card 5 
commences with the application of a time-frequency transform to encode the 
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signal into sub-bands. Essentially, for each track of the acquired, sampled 
and digitized frame, the signal is decomposed into a series of frequency 
bands or sub-bands using a time-frequency transform, such as a wavelet 
transform, and in particular by means of discrete wavelet packet 
5 decomposition. Filtering is performed by the bank of integrated filters 7. The 
number of bands or sub-bands into which the digitized signal is decomposed 
depends on the number of consecutive filtering operations used. In the 
explanation below, reference will be made to a decoding procedure into eight 
bands or sub-bands, which will be indicated below as band 0, band 1, ... band 

10 7. However, as shall be apparent hereunder, this value is not binding; on the 
contrary, an innovative aspect of the method of the present invention consists 
in the capacity to extract from the radio frequency ultrasonic signal a larger 
amount of information for tissue characterization by decoding into a number 
of bands or sub-bands of various width (that is, with amplitude varying from 

15 band to band) and in a number that may be different than eight. 

The filter coefficients can be defined by using a Daubechies 1 6 wavelet. 
By way of example, Figure 6C shows the variation of the Daubechies 16 
wavelet and Figure 6D shows the variation of the DWPT coefficients obtained 
by filtering the digitized and sampled RF signal in Figure 6B. However, the 

20 type of time-frequency transform used may also vary as a function of the type 
of structure or tissue to analyze. For example, different wavelets to the 
Daubechies 16 can be used. Typically, to examine the carotid artery, a 
biorthogonal 3.7 wavelet can be used for characterization of carotid plaques, 
as shall be explained below with reference to some experimental data. 

25 The number of filter coefficients is generally low, for example 32 for the 

Daubechies 16 filter, and the coefficients can therefore be implemented 
easily in a commercially available circuit, which can normally provide filters 
with a limited number of coefficients, not exceeding 64 for each section, with 
decimation by 4, if the filter is not of the linear phase type, while a maximum 

30 of 128 can be obtained with two sections in cascade. The coefficients of the 
filters which provide a wavelet transform are defined recursively and are 
known. For example, it is possible to use the coefficients defined in US-A- 
6066098, the content of which is incorporated in the present description. 

Figure 7A schematically shows a frame acquired by the scanning head 
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of the ultrasound scanner. The frame, referred to below as the "input frame", 
consists of a plurality of n tracks. Each track consists of an analog signal 
represented by the wavefront reflected or back-propagated by the tissues in 
which the ultrasonic excitation signal is propagated. After sampling and 
5 digitization, a frame is obtained (Figure 7B) which also consists of n tracks, 
each of which contains a series of digitized signal samples. This frame will be 
referred to as the "digitized frame". 

By applying the time-frequency transform to the digitized frame, that is, 
by filtering the sampled and digitized signal consisting of the set of 

10 coefficients forming the digitized frame (shown schematically in Figure 7B), a 
distribution in bands or sub-bands is obtained: in practice, for each track from 
1 to n there is, for each instant of time, i.e. for each sample in the time 
direction (that is; in the direction of propagation of the signal), a coefficient in 
each of the bands shown schematically in Figure 7C and indicated as band 0, 

15 band 1 , ... band 7. 

In the example shown the bands are eight and are all of the same 
width. Nonetheless, as mentioned above, it must be understood (as will also 
be explained below with reference to some experimental data) that the 
number of bands into which the signal is decomposed may be greater. With a 

20 greater number of decoding sub-bands from the radio frequency signal a 
greater amount of information on its spectral characteristics, and therefore on 
the structure under investigation, is obtained. Moreover, subdivision into sub- 
bands can also be obtained with frequency intervals (that is, band 
amplitudes) which are not constant, that is, in bands or sub-bands with a 

25 variable amplitude and with filters of any form. The choice of the type of 
decomposition will follow experimental criteria as a function of the type of 
tissue or other structure to be subjected to ultrasound investigation. 

The matrix containing the coefficients of the time-frequency transform is 
indicated below as "transform matrix"; the coefficients of said matrix are 

30 indicated as "spectral coefficients". This is a three-dimensional matrix, the 
dimensions of which are given by: number of tracks, number of samples per 
track, number of bands or sub-bands into which the radio frequency signal 
has been subdivided. 

The subsequent operation consists in extracting the local estimators of 
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the signal for each point of the matrix of the spectral coefficients (and 
therefore for the various points of the input frame from which the matrix of the 
spectral coefficients was obtained with the operations explained above). 

For this purpose, for each temporal instant, that is, for each point along 
5 the dimension of time (that is, of depth) the coefficients of the various bands 
0-7 are taken along each individual track and an interpolating polynomial 
which approximates the variation of these coefficients is determined. This 
operation is shown schematically in Figures 7D, 7E. More specifically, Figure 
7D shows the spectral coefficients in the various sub-bands, relating to track 

10 1 and to track n, by way of example of a generic temporal instant pk. All the 
coefficients of the various bands 0-7 in position f,k and n,k are shown aligned 
with one another respectively in the first and in the second diagram in Figure 
7E, where Ph indicates the interpolating polynomial, that is, the polynomial 
which optimally approximates the spectral coefficients aligned at the depth k 

15 on track 1 and Pl n indicates the interpolating polynomial relating to the 
spectral coefficients aligned at the depth k on track n. 

One or more coefficients of this polynomial constitute a local estimator 
which allows, as will be explained below, extraction of the required 
information on the spectrum of the radio frequency signal. In fact, this 

20 parameter is indicative of the variation of the spectral characteristics of the 
signal. This is because it is related, by definition, to the properties of the 
polynomial which approximates the variation of the spectral coefficients in 
each band for each point of the matrix of the spectral coefficients and 
therefore, indirectly, of the input frame. In other terms, the local spectral 

25 estimator extracted from the interpolating polynomial PI provides a qualitative 
indication of the variability of the spectral characteristics of the RF signal at a 
given point of the input frame. 

Since the interpolating polynomial is calculated at each point of the 
input ultrasound frame, the result of this operation - for each coefficient of 

30 the interpolating polynomial - will be a matrix consisting of the set of spectral 
estimators. This matrix, referred to below as the "local estimator matrix" or 
"local estimator frame", is shown in Figure 7F for a single coefficient of the 
interpolating polynomial. Mere the individual values of the local estimator are 
indicated by ag, where i varies from 1 to n, n being the number of tracks in the 
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matrix, and j varies from 1 to m, m being the number of coefficients per track. 
It must be understood that, depending on the significance of the individual 
coefficients of various orders of the interpolating polynomial, it would be 
possible to create a matrix of local estimators for each coefficient of the 
5 polynomial. Therefore, for a fourth order interpolating polynomial it would be 
possible to construct up to five local estimators matrices, one for each of the 
five coefficients of the polynomial A higher number of matrices can be 

constructed using a polynomial of a higher degree. Not all the coefficients of 

the polynomial may necessarily be useful in the spectra! characterization of 
10 the individual ultrasonic signal. The greater or lesser significance of these 
coefficients is determined experimentally as a function of the type of tissue or 
structure to be investigated. Some non-limiting examples will be described 
below. 

The least squares method of approximation is useful for the calculation 
15 of the values of the coefficients of the polynomial. This method can be used 
to determine the polynomial of degree m < n-1, given a number n of 
experimental observations y 0 , yi, , y n -i corresponding to the abscissae x 0 , 

X1 Xp-1,: 

P(x)-a 0 + ai x+ + a m x m 

20 which gives the best approximation of the detected data. 

Considering the sum S of the squares of the differences between the 
approximate and the observed values, 

H-l 

S = 2> 0 +a, Xj + + a m xj -yjf 

j=o 

the coefficients must be selected in such a way that S is minimal, in order to 
25 make the polynomial approximate the data as closely as possible. By deriving 
S with respect to the coefficients and equalizing the result to zero, and given 

s k = ^T l Xj and v k = 2j x /^j tne s y stem °f equations is obtained, which can 

be written in matrix form as: 

« = £~V 

30 where a is the vector of the coefficients of the polynomial P(x). 

It was observed for the various types of tissue and for various 
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pathological conditions that the values of the coefficients of the interpolating 
polynomial take statistically determined values. By defining the classes of 
values, it was further observed that the statistical distribution, that is, the 
distribution histogram of each coefficient of the interpolating polynomial, has 
5 a variable shape depending on the type of tissue and any pathological 
conditions it may have. According to one aspect of the present invention, the 
method is based on this finding and uses the shape parameters of the 
histogram to extract a greater number of spectral information from the 
ultrasonic signal. 

10 To understand the importance and extent of this variation of the shape 

of the histogram in various conditions and for various tissues, reference 
should be made to Figures 8A to 17B. 

In analyzing the carotid artery to check for the presence and 
consistency of any carotid plaques, for example, it was found that by using a 

15 decoding into sixteen bands with a biorthogonal 3.7 wavelet and by using a 
fourth order interpolating polynomial, not only the third but also the other 
coefficients of the polynomial have a significant distribution for 
characterization of the investigated tissues. It was also found that analysis of 
the distribution histograms show a different morphology of the distribution 

20 histograms depending on the various types of tissue examined. 

Figure 8A shows the B-mode ultrasound image of the carotid artery 
used to obtain the experimental results discussed below. Also identified on 
the image are three windows centered on the lipid plaque, in the free region 
of the vessel, that is, on the blood, and on the calcifications. The coefficients 

25 of the matrix of local estimators were calculated in these regions having 
tissue structures differing from one another, showing that these estimators 
have a statistical distribution which varies according to the type of tissue 
investigated. 

Figures 9A, 9B and 9C represent the variations of the three fourth order 
30 interpolating polynomials relating respectively to plaque, blood and 
calcifications. In other terms, following the digitized radio frequency signal 
decoding and processing method described above, the fourth order 
polynomials which optimally approximate the spectral coefficients (that is, the 
coefficients of the matrix of the transform) have the variation in Figure 9A for 
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the portions of image relating to a lipid plaque, the variation in Figure 9B for 
the portions of image relating to the free region of the vessel, in which blood 
flows, and the variation in Figure 9C for the portions of image relating to the 
regions in which calcification is present. 
5 Figures 10, 1 1 and 12 show the statistical distribution histograms, in the 

various classes of reference, of the five coefficients of the interpolating 
polynomials in Figures 9A, 9B and 9C respectively, obtained by analyzing the 
various regions of a portion of carotid artery. The classes of reference are 
indicated on the abscissas and the population frequencies on the ordinates. 

10 For each of the five coefficients taken into consideration (indicated here 

with AO, A1, A2, A3, A4), three histograms were constructed corresponding to 
three distinct windows positioned respectively in a region of the lipid plaque 
(Figure 10), in a free 'region of the vessel, that is, filled with blood (Figure 11) 
and in a region subject to calcification (Figure 12). It can be seen from the 

15 histograms in Figure 11 that distribution of the coefficients of all degrees is 
very narrow in the free region of the vessel. The relative distribution 
histograms have extremely limited dispersion. On the contrary, in the regions 
affected by lipid plaque and calcification (Figures 10 and 12), the second 
order coefficient (A2) in particular has essentially different statistical 

20 distribution compared with that of the same coefficient in healthy tissue 
(blood, Figure 1 1) and, moreover, essentially different in the two cases of lipid 
plaque and calcification. In other terms, by analyzing the variation of the 
histograms representing the statistical distribution of the second order 
coefficient, it is possible to distinguish whether the region examined contains 

25 a free vessel (blood, Figure 1 1 ), lipid plaque (figure 1 0) or calcification (Figure 
12). 

This means that lipid plaque and calcification introduce different 
spectral modifications in the back-propagated signal, which can be detected 
using qualitative evaluation of the distribution histogram of the second order 
30 coefficients of the interpolating polynomial which approximates the variation 
of the spectral coefficients. 

Analogous considerations can be made regarding the variation of the 
histogram relating to the third order coefficient, which has a different 
distribution in the three windows taken into consideration, centered 
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respectively on lipid plaque, blood and calcification. The histograms relating 
to distribution of the first order coefficient differ to a lesser, although still 
significant, extent. 

Using a decoding into a high number of bands, sixteen in the example 
5 of the data indicated in these figures, use can be made of interpolating 
polynomials of a degree higher than the fourth, for example a fifth order 
interpolating polynomial. From this, even more significant experimental 
results are obtained relative to the statistical distribution of the various 
coefficients. Figures 13A, 13B and 13C show the variations of the fifth order 

10 interpolating polynomials which approximate the variation of the spectral 
coefficients in the case of the carotid artery, obtained with a decoding into 
sixteen sub-bands using a biorthogonal 3.7 wavelet. 

Figures 14, 15 and 16 show the histograms representing the statistical 
distribution of the coefficients (A0-A5) of the various orders of the 

15 interpolating polynomial for the same windows indicated in Figure 8A, 
centered respectively: on the lipid plaque, on the blood and on a region 
affected by calcifications. 

By comparing Figures 9A, 9B, 9C and 13A, 13B, 13C a different 
variation of the interpolating polynomial can be observed especially in the 

20 lipid plaque and calcification region. A different variation of the histogram of 
the coefficients of the interpolating polynomials in the various regions 
corresponds to this, with an even greater difference in statistical distribution 
between the various types of tissue than the difference found in the 
histograms relating to the fourth order polynomial. 

25 Similarly interesting results were obtained in the ultrasound examination 

of the prostate gland. The B-mode ultrasound image used to obtain the 
experimental results set forth and discussed below is shown in Figure 8B. It 
also indicates the window centered on the region of tumoral tissue examined. 
Also in the case of the prostate, using a fourth order interpolating 

30 polynomial and a decoding of the digitized radio frequency signal into sixteen 
bands or sub-bands, a statistical distribution of the first, second, third and 
fourth order coefficients of the interpolating polynomial is noted, which differ 
between healthy tissue and neoplastic tissue, also for the shape of the 
distribution histograms. Figures 17A and 17B show the population frequency 
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histograms of the various classes of reference for the coefficients a1 and a2 
(Figure 17A) and for the coefficients a3 and a4 (Figure 17B) respectively of 
the first, second, third and fourth order, both for the tumoral tissue and for the 
healthy tissue, and moreover in the case of absence of tissue (noise 
5 measured in a sample of water). By comparing the histograms, it can be seen 
that all the coefficients provide information useful for identification of the 
tumoral tissue, as the forms of the histograms of the various coefficients 

change between healthy and neoplastic tissue for all four coefficients. 

These experimental results show on the one hand the relevance (for the 

10 purposes of spectra! analysis) of different coefficients of the interpolating 
polynomial, and on the other the importance of the shape factor of the 
distribution histograms in the various classes of value of the local estimators 
for the purpose of identification of tumoral tissues or in any case of particular 
characteristics in the organs under examination. 

15 Therefore, according to a first aspect, according to the method of the 

present invention shape factors of the distribution histograms of one or more 
of the coefficients of the interpolating polynomial are used to obtain an 
information that can be represented overlaid on the normal B-mode image of 
the ultrasound scanner. 

20 According to a different aspect, the method of the present invention is 

characterized by the (combined) use of several local estimators, that is, of 
several coefficients of various orders of the same interpolating polynomial. 

A possible processing method is described below with reference to 
Figures 7G-7I. The example considers a single matrix of local estimators, that 

25 is, the values of a single coefficient of the interpolating polynomial. The matrix 
of local estimators taken into consideration can be the one relative to any of 
the coefficients of the interpolating coefficient, which was significant in the 
investigation of the type of tissue under examination. The procedure 
described below can be repeated for each of the matrices of local estimators 

30 of interest, each matrix being obtained by collecting the values taken, in the 
ultrasound image, by a specific coefficient of the interpolating polynomial. It 
would also be possible to start with a matrix of local estimators composed of 
a combination of matrices relative to coefficients of a different order of the 
same polynomial. For example, it would be possible to produce a matrix of 
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local estimators which collects the second order coefficients and a different 
matrix of local estimators which collects the third order coefficients and then 
from these a matrix of local estimators in which each coefficient is the sum of 
the homologous coefficients of the first two matrices mentioned. In a different 
5 embodiment, it would be possible to use, as local estimators, the products of 
the local, estimators collected in two or more matrices of local estimators. It 
would also be possible to use all the coefficients of the interpolating 
polynomial. In the case of a fourth order polynomial, it would be possible to 
combine five coefficients (adding or multiplying them, or subjecting them to 

1 0 another type of processing) to obtain a single matrix from the five matrices. 

Returning to the example illustrated, which uses a single coefficient of 
the interpolating coefficient and therefore a single matrix, starting from the 
matrix of the local estimators in Figure 7F, the distribution histogram of the 
values of the local estimators is determined in various classes of reference. 

15 A window with suitably defined dimensions I x J is made to move over 

the matrix of the local estimators taken into consideration. The window is 
positioned in a first position (Figure 7G) and the distribution histogram of the 
coefficients ay of the matrix of local estimators which fall within the window is 
calculated in this position. 

20 As was mentioned with reference to Figure 10 and following, relative to 

the experimental results, the distribution histograms of the coefficients of 
various orders of a given polynomial vary to take different shapes depending 
on the type of tissue found in the specific portion of image examined, 
enclosed by the window. Once the investigation window with dimensions IxJ 

25 has been positioned in a, given position on the matrix of the local estimators, 
the distribution histogram is performed and a shape factor, for example the 
standard deviation a, of the histogram is calculated. By moving the window 
one step at a time both along the rows and along the columns of the matrix of 
local estimators, a series of standard deviation values equal to the number of 

30 coefficients of the matrix of local estimators is obtained. By indicating with try 
the value of the standard deviation relating to the distribution histogram of the 
coefficients of the matrix of local estimators calculated with the window in the 
position ij, it is possible to construct a matrix with a number of coefficients 
equal to the number of coefficients ay of the matrix of local estimators, but the 
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value of which has been modified and processed taking into consideration 
the shape parameter (in the specific example the "standard deviation erg) of 
the distribution histogram. This matrix will be composed of m x n values by 
bjj = f(a]j, ag) 

5 where f is a generic function of ay and a|. For example, each coefficient ay of 
the matrix of local estimators could by divided by the corresponding value of 
CTij obtaining: 

- by g a i /g j j 

It must, however, be understood that the type of function that relates 
10 the coefficients by to the coefficients ay and to the shape factor 0% of the 
corresponding distribution histogram may differ from the one indicated. For 
example, it would be possible to perform a multiplication rather than a 
division, or a different combination of these values. 

Moreover, it would also be possible to consider for a plurality of 
15 coefficients of different orders of the interpolating polynomial a corresponding 
plurality of matrices of local estimators and to combine these together, 
weighting them with the shape parameters of the relative distribution 
histograms. 

The shape factor or shape parameter of the histogram which is used 
20 can also differ from the standard deviation. For example, it would be possible 
to use the kurtosis, the symmetry index, higher degree moments or a 
parameter indicative of the most populated classes. In general, the 
coefficients of the matrix of the local estimators could be processed with any 
shape parameter characterizing the distribution histogram of the coefficients 
25 of said matrix in a specific window, for example centered on the position of 
each individual coefficient of the matrix of local estimators. 

The matrix of dimension m x n which collects the coefficients by is 
hereunder called matrix of weighted local estimators. This means that each 
coefficient of said matrix has been "weighted" by a shape factor of the 
30 distribution histogram of the corresponding window, where "weighted" is 
intended as the fact that this coefficient has been combined in some way with 
a shape parameter, such as standard deviation or another suitable 
parameter. 
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It must be understood that the operation described may be performed 
separately for all or for some of the coefficients of the various orders of the 
interpolating polynomial which are useful for spectral characterization of the 
ultrasonic signal. The operation performed on the individual matrix of the 
5 local estimators which collects the coefficients of a generic order of the 
interpolating polynomial can be the same for all the orders taken into 
consideration, although it would also be possible to process the coefficients 
of the various orders of interpolating polynomial in different ways. For 
example, the third order coefficients could be weighted by dividing by the 

10 parameter of standard deviation, while the second order coefficients could 
either be weighted, or another operation of combination with a different 
parameter, such as kurtosis, or an operation differing from dividing the local 
estimator and the shape parameter could be performed. 

In any case, the general concept is to combine the value of the 

15 coefficients representing the local estimator with a coefficient that takes into 
account the shape of the distribution histogram of said coefficients in a 
specific window. 

Figure 7J schematically indicates the procedure which can be 
performed when, instead of a single coefficient of the interpolating 

20 polynomial, all the coefficients of said polynomial are used to obtain a more 
complex and elaborate set of weighted local estimator. The matrix of the local 
estimators is therefore a three-dimensional matrix with dimensions m * n * k, 
where m is the number of samples per track, n is the number of tracks in the 
digitized and sampled ultrasound frame and k is the number of coefficients of 

25 the interpolating polynomial (polynomial of degree k-1). The coefficients of 
the matrix are indicated with ag (1) , ... ai/ k) . 

The coefficients af>, ...a 0 (k) are combined with one another to 
obtain a coefficient cy of a composed local estimators matrix. This matrix is a 
two-dimensional matrix with dimensions n*m. Each of its coefficients is given 

30 by 

c^f(a g < 1 > a;/ 2 >, ...aft 
with f generic function. For example, it may be a summation or a product of 
all the k coefficients of position ij of the three-dimensional matrix. 

On this matrix eg the shape coefficients ay, for example the standard 
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deviation, of the distribution histograms are calculated by means of the 
procedure described above with reference to Figures 7G, 7H. The matrix of 
the weighted local estimators, again indicated with by, and from this the 
population matrix % are obtained from the shape coefficients cy and from the 
5 shape parameters ct s . 

The procedure is similar in concept to the one described with reference 
~to Figures 7D-7H, although with the difference that coefficients (weighted 
local estimators) which contain information deriving from all the coefficients of 
the interpolating polynomial, instead of only one of these, are contained in the 
10 final matrix. 

Intermediate solutions, in which only some of the coefficients of the 
interpolating polynomial are used, are also possible. 

In a variant of embodiment, it would be possible to first perform 
weighting of the individual coefficients ay through the shape coefficient ay and 
15 only subsequently combine the weighted local estimators with one another. 
Essentially, in this case after determining the matrix of the spectral 
coefficients the interpolating polynomials are determined and from their 
coefficients the matrix of the estimators ay (k> . For each value of k the 
statistical distribution, that is, the shape of the histograms, is determined and 
20 a shape factor, for example the standard deviation a-^, is calculated. The 
weighted local estimators are then calculated 

by< k ^ay< k U/ k) ) 

Finally, the weighted local estimators are combined to obtain 

Cy = f(bi/ 1) by (k) ) 

25 where f and f are generic functions. 

According to a different aspect of the invention, the procedure does not 
weight the coefficients or local estimators. This different embodiment of the 
procedure does not allow the use of information contained in the shape of the 
distribution histograms, but still allows the use of the information content of 

30 several coefficients of the interpolating polynomial. Figure 7K schematically 
represents the processing procedure in this case. Identical symbols indicate 
identical or equivalent elements to those in Figure 7J. The matrix cy is no 
longer processed to obtain the shape factor cry, but its coefficients (local 



-25- 



rei/ITZUU4/UUU54S 



estimators, not weighted) are subjected to statistical analysis to directly 
obtain the population matrix, indicated here with fy'. 

Having obtained the matrix of the weighted local estimators, the values 
of the coefficients of this matrix are classified according to classes of values 
5 which have been found, experimentally, to be related bi-univocally to a 
specific tissue structure. In other terms, it is experimentally proven that a 
certain type of tissue (for example the neoplastic tissue of the prostate) 
modifies the ultrasonic excitation signal so that a certain weighted local 
estimator takes values that fall within a certain class of values. The 

10 ultrasound scanner is then programmed to perform, on an input frame, the 
operations described above to obtain the matrix of weighted local estimators. 
The coefficients of the matrix thus obtained are then compared with the class 
of values characteristic of the neoplastic tissue and a color image formed of 
points of the matrix of the weighted local estimators, the values of which fall 

15 within the characteristic class of the neoplastic tissue, is overlaid on the black 
and white B-mode image. This color indication provides the ultrasound 
operator with the required information, showing the presence and position of 
the neoplastic tissue. 

It is understood that classification can be performed for various types of 

20 tissue or structure. For example, if it is desired to show the presence of 
healthy tissue, it is enough to identify the class of values of a certain weighted 
local estimator which is characteristic of the presence of healthy tissue. 
Likewise, when it is desired to simultaneously or sequentially show on the 
same image the presence of two different types of tissue or structure (for 

25 example lipid plaque and calcification in the carotid artery), it is enough to 
identify the two classes of values which characterize these two types of 
structure and then filter the coefficients of the matrix of the local estimators 
with the characteristic values of these classes and present the result 
simultaneously or sequentially on the B-mode ultrasound image, if necessary 

30 and preferably with two different colors. 

To experimentally identify the classes of values which for a certain local 
estimator characterize a specific type of tissue, it is possible to proceed by 
examining an experimental sample in which the structures of interest are 
identified by means of histological analysis. In the case of the prostate, for 
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example, the homogeneous regions of healthy tissue, adenocarcinoma and 
adenofibroma are identified. The radio frequency ultrasonic signal, obtained 
in response from scanning of the experimental sample, as described above, 
is then processed to obtain a single matrix of weighted local estimators, 
5 Due to the way in which the matrix of local estimators has been 

obtained, a point of the B-mode image corresponds to each coefficient of the 
matrix. Therefore, as the region is known, in which a certain structure (for 
example the adenocarcinoma) is localized, statistical distribution of the local 
estimators is determined En these regions and the class with the greatest 

10 population frequency is identified. That class will be associated, for that 
specific organ, with the presence of adenocarcinoma. Therefore, each time 
an ultrasound analysis is performed on an analogous organ, the weighted 
locai estimators will be classified and those falling within the class of 
reference of the adenocarcinoma will be highlighted by overlaying a color 

15 scale on the B-mode image in which the pixels, corresponding to the regions 
of the image in which the values of the weighted local estimators fall within 
the characteristic class of the adenocarcinoma, are colored. An analogous 
procedure will be used for the other types of tissue (healthy tissue, 
adenofibroma, in the case of the prostate). 

20 More detailed processing envisages working on at least one pair of 

matrices of local estimators both in the phase to identify the classes of 
reference, and in the phase to analyze an organ subjected to ultrasound 
investigation. One example of processing of this type, based on three distinct 
local estimators, is now described with the aid of figures 23A and 23B. 

25 Acquisition of an ultrasound frame is once again performed on the 
experimental sample, which is then processed as described above to obtain 
the three matrices of weighted local estimators. 

A first two-dimensional graph shows, respectively on the abscissas and 
on the ordinates, all the values of the coefficients of two matrices of weighted 

30 local estimators, for example those obtained by order 0 and order 1 
coefficients of the interpolating polynomial. A second graph shows the 
weighted local estimators obtained by the order 0 and order 2 coefficients. 
Figure 23A shows, on the abscissas, the coefficients of the matrix of 
weighted local estimators obtained by the first order coefficients of the 
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interpolating polynomial for a prostate. The coefficients of the matrix of 
weighted local estimators obtained by order 0 coefficients of the same 
polynomial are shown on the ordinates. Figure 23B shows, on the abscissas, 
the weighted local estimators obtained by second order coefficients and, on 
5 the ordinates, those obtained by order 0 coefficients. 

Therefore, a pair of values (the coordinates of the point), composed by 
the value of one and of the other of the weighted local estimators used, is 
associated with each point of each of the two graphs. 

The graphs obtained in this way can be overlaid on the input ultrasound 
10 frame. Therefore, it is possible to associate with each point of the input 
ultrasound frame pairs of classes of values of the weighted local estimators 
examined and corresponding to this point. 

By means of histological analysis, the regions of the tissue which are 
homogeneous in nature (healthy, neoplastic, etc.) are identified on the 
15 ultrasound frame. In Figures 23A and 23B these regions are distinguished by 
three different shades of grey. As a pair of values of the weighted local 
estimators correspond, in both graphs, to each point of the ultrasound frame, 
it is possible to associate a class of values for said weighted local estimators, 
with each homogeneous region of the tissue examined, said class being 
20 univocally associated with that specific type of tissue. 

In order to identify each region bi-univocally, it is important to select a 
set composed of pairs or triplets of classes of values of the weighted local 
estimators, that are not overlaid on any other associated with a different 
region. The operation to extract these sets can be performed through, for 
25 example, OR and AND logic operations, 

As mentioned, the population frequency of the sets thus identified is 
represented with a color scale, different for each set, on the conventional B- 
mode image and is directly proportional to the intensity of the population 
frequency. 

30 Experimentally, as can also be noted in Figures 17A, 17B, the regions 

which is desired to show (in the example the neoplastic tissue) can normally 
be highlighted by combining simply with a sum the images obtained from all 
the coefficients of the interpolating polynomial, which can be highlighted by 
combining, in this case, with an AND operation the sets of values obtained 
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from all the coefficients of the interpolating polynomial. 

The procedure described above for experimental correlation between 
tissue characteristics detected by histological analysis and classes of values 
for the weighted local estimators is also valid in the case of non-weighted 
5 local estimators, for example obtained by combining different degree 
coefficients of the interpolating polynomials with one another. 

In the experimental results discussed above and indicated in Figures 
9A-17B, decomposition of the RF signal into sixteen bands or sub-bands of 
the same width was used. This decomposition provided useful results both in 
10 examination of the prostate, and in examination of the carotid artery. These 
two types of organ were chosen as they have tissue structures with different 
dimensions. The prostate gland typically has structures with dimensions 
larger than the carotid artery. The different dimension of the structures that 
reflect the ultrasonic excitation signal cause different distribution of the 
15 frequencies in the return signal. Subdivision into sixteen bands of the same 
amplitude makes it possible to obtain a sufficiently significant resolution of the 
RF signal returned from both the organs examined. 

Nonetheless, in general it will be understood from the above description 
that the choice of decomposition into bands (number and amplitude of bands) 
20 may differ according to the type of tissue or structure to be analyzed. The 
smaller the structures to be studied in the analyzed organ are, the greater the 
resolution of decomposition into high frequencies will be, and vice versa. This 
is confirmed by the experimental results discussed below with reference to 
Figures 18 to 21, obtained again from processing an ultrasound image of a 
25 portion of carotid artery. 

Figures 19, 20 and 21 show the B-mode image of the portion of carotid 
artery analyzed. The digitized and sampled RF signal was decoded using a 
discrete wavelet into six bands of variable amplitude. Figure 18A shows the 
tree representing the filter process used and Figure 18B shows the six 
30 decomposition bands. It can be observed that the amplitude of the bands 
increases from low frequencies to high frequencies. The sampling frequency 
of the signal is indicated with fc. Figure 18C shows the variation of the 
interpolating polynomial obtained applying the aforesaid decomposition to the 
B-mode image in Figures 19, 20 and 21. 
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Figure 19 shows, in the top part, the B-mode image of the carotid artery 
examined, below this the variation of the fourth order interpolating polynomial 
(corresponding to Figure 18C) and subsequently the distribution histograms 
of the coefficients A0-A4 of the various degrees of interpolating polynomial, 
5 calculated at the level of a window (shown overlaid on the B-mode image) 
centered on the lipid plaque. Figure 20, with the same layout of the images, 
shows the interpolating polynomial and the distribution histograms obtained in 
a window centered on the free region of the vessel, that is, on the blood. Also 
in this case the window is shown overlaid on the B-mode image of the carotid 
10 artery. Finally, Figure 21 shows, with the same disposition, the result obtained 
in a window centered on a portion of the B-mode image in which calcification 
phenomena are recognized. 

By comparing Figures 19, 20 and 21 it can be seen that with 
decomposition into six bands of amplitude increasing from the low to the high 
15 frequencies, the variation of the interpolating polynomial does not change 
appreciably in the three regions observed (lipid plaque, blood and 
calcifications). This is a first indication of the fact that the decomposition used 
does not provide resolution appropriate for examination of this type of organ. 
This is due to the fact that the dimensions of the structures to be recognized 
20 are small and, therefore, the spectral information they supply is located at the 
high frequencies of the response spectrum, that is, where the decomposition 
used requires bands with a greater amplitude. Contrary to this, the results 
obtained on the same organ with decomposition into sixteen bands of the 
same amplitude (Figures 9A to 16) show that the resolution obtained with this 
25 decomposition is sufficient to provide significant results. 

From the distribution histograms in Figures 19, 20 and 21 it can be 
noted that in fact (although the variation of the interpolating polynomial is 
essentially the same in the three regions of lipid plaque, blood and 
calcification) the shape of the distribution histograms of the various orders of 
30 coefficients changes. Nonetheless, the difference between distribution 
histograms of homologous coefficients in the three investigated regions is not 
as marked and significant as the difference obtained using signal processing 
based on decoding into sixteen bands of identical amplitude. In conclusion, 
comparison of the results obtained with the two types of decoding into sub- 
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bands shows that the choice of suitable decoding, with regard to number of 
bands and their amplitude, allows the method of analysis of the present 
invention to be adapted to obtain optimal results according to the type of 
tissue or structure investigated. The type of decoding or decomposition can 
5 be selected experimentally, as can the most suitable degree of the 
interpolating polynomial, its most significant coefficients and the most 
relevant shape coefficients of the distribution histograms of said coefficients. 

Figure 22 shows as an example a decomposition tree of the spectrum 
of the radio frequency signal by the Discrete Wavelet Packet Transform 
10 (DWPT). According to the type of ultrasound investigation to be performed 
and to the characteristics of the examined structure, the investigation method 
can use any one of the various decomposition possibilities along the tree in 
Figure 22. 

Clearly, the drawing shows only one practical embodiment of the 
15 invention, which can be varied in its forms and arrangements without 
departure from the scope of the principle of the invention. The presence of 
any reference numbers in the attached claims has the sole purpose of 
facilitating reading with reference to the preceding text and the attached 
drawings, and does not limit the scope of protection of the claims. 



